ON CONVERGENCE IN SMOOTH GRADIENT 
SYSTEMS WITH BRANCHING OF EQUILIBRIA 



V.A. GALAKTIONOV, S.I. POHOZAEV, AND A.E. SHISHKOV 

Abstract. Our basic model is a semilinear elliptic equation with a coercive nonlin- 
earity, A-0 + f{ip) = in S7, = on 917, where il C is a bounded smooth domain. 
Our main hypothesis (Hr) on the resonance branching is as follows: if a branching of 
equilibria occurs at any point ip with a /c-dimensional kernel of the linearized operator 
A + f'{ip)I, the branching subset 5*^ at ■0 is a locally smooth fc-dimensional manifold. 
Using the corresponding parabolic flow 

ut = Aw + f{u), 

with bounded initial data wq, we then prove that, under the (Hr), the subset of equilibria 
$ = {ij)} is evolutionary complete, i.e., any evolution orbit converges as t ^ cxi to a 
single element of $. We also treat the non-coercive case f{u) = \uY'~^u for p g (1, jfZ^), 
where the elliptic problem is known to admit at least a countable subset of solutions {V'fc} 
due to the classical Lusternik-Schnirel'man category theory. The results are extended 
to higher-order elliptic operators — (— A)™!* -I- lul^*"^!* with Dirichlet conditions on dVt, 
1 < p < TO > 2, with the corresponding adaptation of the parabolic flows. 

For = 1, the first result on stabilization to a single equilibrium is due to T.I. Ze- 
lenyak (1968). We show that Zelenyak's approach based on Lyapunov's function analysis 
can be extended to general gradient systems in Hilbert spaces with a similar smooth res- 
onance branching. We also cover the case of their asymptotically small non-autonomous 
perturbations. In general, the developed approach represents an alternative (and im- 
proved) method on stabilization in Hale-Raugel (1992) and other later similar techniques 
in gradient system theory. 



Dedicated to the memory of Professor T.I. Zelenyak 



This is an extended version of the paper |T2], containing extra comments and updated 
references. 

1. Introduction: the convergence problem and evolution completeness 

1.1. Semilinear elliptic equations and the convergence problem. Let Q C 

be a bounded domain with the smooth boundary dQ. Our basic model is the classical 
semilinear coercive elliptic equation from combustion theory 

(1.1) V'ii/j) =A^p + fiij) = in fi, ^^J = on dil, 
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where / is a given function satisfying 

(1.2) / e C\M.), f{u) < for \u\ > 1. 

This non-optimal coercivity condition is sufficient for demonstrating the approach. For 
instance, f{u) = u — (an analytic nonlinearity) or f{u) = \u\^~^u — with any 

q > p > 1. The operator in (11.11) is potential and coercive, so the classical variational 
theory establishes existence of solutions as critical points of the corresponding functional 

(1.3) V{ij) = -I J\Dij\^dx + J F{4j)dx, 

where F' = /; see e.g., [20]. If / does not satisfy the coercivity condition (II. 2p . then, in 
the most well-known power case 

(1.4) f(u) = \u\P-^u with I <p <Ps = 

versions of the classical Lusternik-Schnirel'man (L-S) theory [22]-[23] (see e.g., [6] and [281 
[29] ) establish existence of at least a countable subset of critical points of the corresponding 
functional (II. 3p . 

By $ = {-0} we denote the subset of all bounded solutions of the elliptic problem (11.11) 
in both coercive and non-coercive cases. The evolution completeness of the subset $ of 
solutions of (11.11) is defined by introducing the corresponding parabolic equation 

(1.5) ut = V'{u) =Au + f{u) in QxR+, u = on dQ x M+, 

with arbitrary bounded initial data uo. It follows from the classical parabolic theory 
[To] , that under hypotheses (II. 2p the orbit {u{-,t),t > 0} is uniformly bounded. For 
such smooth gradient systems, one can define its tu-limit set uj^uq), which is non-empty, 
compact, invariant and connected in the topology of L^(fi) (or C{Q)), and 

(1.6) uj{uo) C $; 

see [H] and [211 pp. 483-487]. The notion of the evolution completeness of the subset $ 
of nonlinear eigenf unctions was introduced in [llj for the porous medium equation 

Ur = A(|m|™-1m) + (m > 1) 

in a bounded domain. In this case, is not included into $ (as usual, is not an 
eigenf unction) and this determines some specific difficulties of the asymptotic analysis. 

In the present problem, G so that the evolution completeness of $ just means that, 
for the parabolic flow (II. 5p . 

(1.7) for any data Uq, 3 a -0 G $ such that uj^uq) = {ip}. 

Indeed, this is the classical convergence problem in the theory of dynamical systems. 
Obviously, (11.71) is true for any smooth gradient system, if $ consists of isolated points (a 
straightforward consequence of the connectedness of uj{uo)), i.e., 

(1.8) $ is discrete =^ evolution completeness (II. 7p holds. 

On branching situations: two approaches. The case, where $ is not discrete and 
consists of continuous families due to branching of equilibrium points was for a long time 
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a problem of concern. Existence of various branches of stationary solutions is established 
by classical variational techniques; see [20] and [19J for applications to semilinear elliptic 
problems. In some cases, the fibering method ^2Bl [29] shows a clear picture of the con- 
nection of the number of branches with the algebraic non-monotonicity of f{u), so, in 
general, $ may contain continuous sub-families; see examples in Section [31 

The first approach to convergence in parabolic problems in the presence of branch- 
ing was formulated in 1968 by T.I. Zelenyak [3l], who proved (11. 7p for arbitrary one- 
dimensional uniformly parabolic equations with smooth coefficients 

(1.9) ut = a{x,u,Ux)uxx + b{x,u,u^) 

with arbitrary nonlinear boundary conditions. Ten years later, in 1978, a similar result 
was proved by H. Matano [25] by Sturm's zero set argument, which turned out to be 
applied to arbitrary smooth fully nonlinear uniformly parabolic equations for = 1 

(1.10) ut = F{x,u,Ux,Ux^). 

Matano's approach is principally one-dimensional, since Sturm's Theorem on zero sets 
is not available in higher dimensions, cannot be extended to parabolic PDEs in for 



Another direction of the development of the stability theory of gradient systems is 
associated with the Lojasiewicz-Simon inequality and technique. It is well-known that, 
for gradient dynamical systems in a real Banach space E 



where V is a functional, with the Frechet derivative V : E ^ E', E' being the dual 
space, the Lojasiewicz-Simon inequality is an effective tool of studying of the stabilization 
phenomena. In particular, it completely settles the case of analytic nonlinearities; see 
references in [5]. Given an equilibrium point il), with V'{il)) = 0, this classical inequality 
has the form 



exponent; see basic references and historical comments in If (11.12!) with a fixed ex- 
ponent 9 holds for the critical manifold ^ Thm. 3.10], then, under natural assumptions 
on the spectral properties of linearized operators (all those hold for our elliptic case), it 
is possible to establish: 

(i) that the cu-limit set ujIuq) of any bounded orbit consists of a unique equilibrium, i.e., 
flTTTll holds, and 

(ii) the rate of convergence to "0 as t ^ oo. 

For instance, in the case of finite regularity, the best possible constant 6 = \ guarantees 
the exponential convergence, [T6l Thm. 1.1], 



> 2. 



(1.11) 



ut = V'{u) for t > 0, u{0) = uoe E, 




(1.13) 



u{t) - i/jWe < ce 



-St 



as t — i> oo (5 > 0). 
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In general, roughly speaking (see examples in [5]) 

(1.14) de{0,l) =^ \\u{t) -tpllE < ast^oo. 

Concerning applications of the Lojasiewicz-Simon approach to parabolic equations like 
(11. 5p . it does not assume any restrictions of the dimension A^, though, in presence of 
branching, checking the validity of the corresponding inequality on the singular subset 
becomes an extremely difficult problem; see examples in [5]. On the other hand, in 
Section [5] we show that the Lojasiewicz-Simon inequality fll.l2p does not apply to general 
C°° dynamical systems, i.e., the exponent 6 G (0, |] does not always exist. Actually, 
this means that this kind of analysis needs other non-rational moduli of continuity in the 
inequality fll.l2p . 

Counterexamples and the resonance branching hypothesis. It is well-known for 
a long time (see [27j) that, for general non-analytic gradient systems on the plane M^, the 
convergence result (11. 7p is not true. The idea of such a construction has been extended 
to prove existence of non- convergent orbits for special parabolic equations like (II. 5p with 
C°° function / = f[u,x)] see [20] and earlier references therein. The principal part of the 
construction in [271 ED] uses the fact that the dimension of manifold of equilibrium points 
is less than the dimension of the corresponding kernel of the linearized operator meaning 
a certain "defect" of dimensions. 

It is well-known that there are special (often similar to ours) cases of branching of 
equilibria, where Hale-Raugel's approach applies to guarantee convergence of the orbits of 
gradient systems; see [l5j, where a survey and further references are given. This approach 
essentially relies on spectral properties of the linearized operator (main Hale-Raugel's 
hypothesis is that has multiplicity at most one, or /c > 1 under special hypothesis) 
and uses properties of stable, unstable and center manifolds, that are difficult to justify 
for some less smooth equations. We also refer to more recent paper [4J, where a similar 
approach to stabilization is used and other references can be found. 

Our method is different and uses Zelenyak's ideas (1968) from one dimension that are 
mainly connected with Lyapunov functions only. It is important that Zelenyak's approach 
yields the exponential rate of convergence to a single equilibrium and can be extended 
to non-autonomous perturbations of gradient systems, for which the classic stable-centre 
manifold theory fails. 

Let us present the main "resonance" hypothesis, under which non-convergent orbits do 
not exist. We formulate it in a general form, where, for the PDEs (11. 5p . by the Frechet 
derivative we mean V"{ijj) = A + /'(t/^)/ assumed to be a Fredholm operator of index zero 
at any equilibrium point from $. 

Hypothesis (Hr). If branching occurs at any equilibrium and dimker ^"(-0) = k > 
1, then the corresponding branching subset Sk of equilibrium points is a locally smooth k- 
dimensional manifold {a relatively open surface). 

As a first comment, we note that, obviously, (Hr) always holds in the linear case with 
f{ip) = Xip, where branching of equilibria occurs for A G o"(A). Then the subsets coincide, 

(1.15) Sk = ker y"(^) (V"(^) = A + A/). 
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Indeed in the linear case the convergence of all the orbits follows from the completeness 
and closure of the eigenf unctions of A as a straightforward consequence of general theory 
of self-adjoint operators, [3], and we do not treat this case here. 

The condition in (Hr,) cannot be improved in the sense that, in a non- resonance branch- 
ing, where the smooth stationary manifold Sk has dimension k < dimker this 
makes it possible to apply the idea of construction of non-converging orbits from [271 [30] . 
According to our proof, any defect of dimensions destroys the approach. Concerning the 
resonance branching condition, we refer to the classical branching theory for nonlinear 
equations in Banach spaces [33l Ch. 7,8], [21 Ch. 6], [H § 30.1] for the conditions, which 
guarantee necessary branching in variational setting, and also to [12] , where applications 
to elliptic problems of interest are discussed in detail. Notice that, for fll.ip . the branch- 
ing parameter is not available explicitly in the equation and such branching cases are 
less studied in the literature. We will present simple examples showing that in standard 
nonlinear elliptic problems, the resonance branching with k = N — 1 actually exists in 
dimensions N > 2. 

1.2. Main result, plan and extensions. We now state the main result of the paper, 
which, as we have mentioned, are established by using key Zelenyak's ideas from \3^. 

Theorem 1.1. Let f be a function. Let $ = {tp} be the subset of all bounded solutions 
of the elliptic problem i\Ll\i and let (Hr) hold. Let, given bounded initial data uq, u{x,t) 
be the unique classical solution of (I i . 5|) . Then: 

(i) in the coercive case there exists a unique ip E ^ such that u{uo) = {if^}; 

(ii) in the non-coercive case (1.4), the same is true for any uniformly bounded orbit; 



(iii) in both cases, the convergence to tp is exponential. 

The hypothesis of boundedness of the orbit in (ii) is essential since (11. 5p admits solu- 
tions, which blow up in finite time and hence exhibit entirely different asymptotic patterns 
(specific asymptotic techniques are explained in [I3|)- As we have mentioned, for analytic 
nonlinearities f{u) (for (11.41) this means p = 3,5, ...), the evolution completeness conclu- 
sion follows from the classical Lojasiewicz-Simon inequality; see references in [311], where 
versions of such inequalities were proposed for some non-analytic settings. The result (iii) 
suggests that in the resonance branching, the Lojasiewicz-Simon inequality with = ^ 
holds in a neighbourhood of 5"^, though we do not study this. 

The layout of the paper is as follows. Theorem 11.11 is proved in Section [2l In Section 
[3] we present some examples showing the possibility of resonance branching. We show in 
Section [5] that, for some C°° nonlinearities, the Lojasiewicz-Simon inequality fails for any 
arbitrarily small 6' > 0, so the stabilization technique needs other non-rational moduli of 
continuity in (11.121) . which are proposed. It follows from (iii) in Theorem 11.11 that such 
essentially non-analytic counterexamples correspond to isolated equilibria. 

The results are true for smooth higher-order equations. E.g., Theorem 11.11 is valid for 
a general 2mth-order parabolic flow for any m > 1 

(1.16) ut = {-l)"'+^A'^u + f{u) mQxR+, 
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with m zero Dirichlet boundary conditions. The critical Sobolev exponent for (II. 4p is 
(1-17) Ps = w^ (iV>2m). 

We can also consider quasilinear fourth-order uniformly parabolic equations like 

(1.18) Ut = -Ag{Au) + f{u), g'{s)>co>0, 

or with others potential operators in various metrics. In Section [6] we extend the approach 
to smooth gradient systems in Hilbert spaces. In Section [7] we treat the case of asymp- 
totically small non- autonomous perturbations of gradient systems. We show that the 
main result of convergence (11.71) of bounded orbits remains true for perturbed parabolic 
equations (11. II) such as 

(1.19) Ut = Au + f{u) + h{t)g{u), where h{t) — * as t — > oo, 

and g{u) is a given smooth function, and also applies for similar perturbations of PDEs 
(11.161) and (ll.lSp . We then impose a certain restriction on the rate of decay of h{t) as t — 
CX3, but, anyway, we can treat non-exponentially small perturbations, e.g., h{t) = 0{^). 

2. Proof of Theorem 11.11 

Thus we consider the uniformly elliptic equation (11.11) and the corresponding uniformly 
parabolic flow (11.51) . Notice that the subset of all stationary solutions $ = includes 
the trivial one ip = 0. Obviously, is exponentially stable for ( 11. Sp with nonlinearity ( 11. 4p 
and hence must be taken into account in the evolution completeness analysis. 

Let us mention again that according to the potential structure of the equation (11.10 
and to (Hr), we are going to use the following properties: 

(i) Sk is open, so any ip E has a neighbourhood relatively open in S'^, and 

(ii) Viip) = B = const., for all ip G Sk- 

2.1. Proof Theorem II. li Thus we consider a smooth flow (11.50 in, say, Hq{Q) and 
study the a;-limit set uj{uq) of a uniformly bounded orbit {u(t), t > 0}. Take an arbitrary 
equilibrium i/j from the w-limit set (jj{uo), which is non-empty, compact and connected for 
such smooth gradient systems, So that there exists a monotone sequence {tk} oo 
such that 

(2.1) u{tk) ip as A; ^ CX3. 

First of all, if ^ a(y"(^)), where V"{ij) = A + f'iij)I, then ip is isolated [331 |2Q] and 
the convergence result (11.70 follows for given uq. Second, assume that, under hypothesis 
(Hr), G a{V"{ip)) and ip is not isolated and belongs to a smooth manifold Sk branching 
at i/j, where k coincides with the dimension of ker V"{iIj). 

Our further analysis uses the structure of Zelenyak's proof in the case = 1, [M] and 
consists of seven steps. 

Step 1: optimal approximation of u{t) on Sk- Let Sk = {"0^}, where G denotes 
local coordinates on Sk such that i/j = ipo. Given a point u{t) for some t S> 1, when 
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u{t) ~ ip, we define the optimal approximation ilj^{t) of u{t) on 5"^ by minimizing the 
distance 

(2.2) hifi) = \\u{t)-ij^\\, 

where || ■ || denotes the L^(f2)-norm induced by the scalar product (■, ■). Then inf^/i(/z) is 
attained at some /i = ~ 0, and, by construction, the orthogonality condition holds, 

(2.3) u{t)-t/j^^t) ±kerV"(^^(t)). 

In the further evolution analysis we deal with equilibria il)^ in 5"^ for /i ?a 0. It suffices to 
prove that there exists the limit 

n{t) — > as t — > oo. 

Step 2: a priori bound for the linearized stationary equation. We linearize the 
stationary equation (II. ip at i/j = ipnit) consider the corresponding inhomogeneous 
problem 

(2.4) V"{ij,^,))w = g, 

where G is a given function. Then, in view of the orthogonality to the kernel, fl2.3p . 
by the standard theory of elliptic self-adjoint operators (see e.g. [2], [3]), there exists a 
constant Ci > such that 

(2.5) \\w\\h2 = II^W - ^mwII^^" ^ C'lllfi'll- 

Step 3: spectral gap. By hypothesis (Hr), G a(y"{tpfj)) for all // ~ 0. Moreover, the 
resonance condition also assumes that the spectral gap 

(2.6) A(/i) = inf{|A| : A G a(A'(V^^)), A ^ 0} 

is uniformly bounded away from zero, i.e., there exists a constant C2 > such that 

(2.7) A(/i) > C2 > for all ^ ^ 0. 

This follows from the fact that, by (Hr), the kernel ker A'('0^) changes continuously with 
/i ~ as the tangent space of the changing continuously points tp^ on the given smooth 
manifold Sk- 

Step 4: estimate for the linearized parabolic equation. We now set w{t) = 
u{t) — ip^{t)- Consider next the parabolic equation (11.50 . linearize the right-hand side 
at M = 'U'^(t) and write down it as the inhomogeneous elliptic equation (12.40 . Then, by 
linearization, we obtain an extra quadratic term gi{u) = | /"(^^)(m — ■?/'^(j))^, so that, in 

m, 

(2.8) g = ut + gi{u), with 1 5(1 (m)| < CgUti (t) - ?/;^(t)| 1^2- 

By convergence in such a smooth parabolic equation, we have that the smallness of 
w{x, t) = u{x, t) — w^(t)(x) implies the smallness of the derivatives, so that we may assume 
that C3 is uniformly bounded on any intervals [tk,tk + T] as k 00. 
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Combining the elliptic estimate (12.51) and the parabolic one (12.81) yields the following 
bound: 

(2.9) \Ht)-ij^^t)\\l2<C4utf, 
for all t ^ 1 such that the orthogonality (12.31) holds. 

Step 5: local exponential convergence of the Lyapunov function. The Lyapunov 
function V{u{t)) is strictly monotone increasing on evolution orbits, 

(2.10) ^^V{u{t)) = J{u,y>0 {u,^0), 

n 

so that there exists the finite limit 

(2.11) V{u{t)) ^5+ as t ^ oo, 

and then V = B on Sk- Fixing the parameter of the optimal approximation /i = /u(t), by 
using the standard expansion at ■0^, one can see that 

(2.12) \V{u) - V{^^)\ < C^Wu - ijJl, < C4u - ijjl,. 
Therefore from fl230D . by (EH) and (127721) . we have that 

(2.13) A [5 _ v{um = - J{u,)' < -C4u - ^,(,)||^. < -Cj[B - V{um. 

n 

Integrating this inequality yields the local exponential convergence of V{u{t)) to B on 
any arbitrarily large bounded intervals t E [tk,tk + T] as A; — > cxd, 

(2.14) B - V{u{t)) <[B- V{u{tk))]e~^'^'~'^\ 

Step 6: local exponential convergence of V{u{t)) implies exponential conver- 
gence of u{t) to equilibrium. Here we use Lemma 4 in [31], actually establishing a 
weighted Gronwall's-type inequality by using a discrete partition technique. Later on, 
finite partitions have been widely used for more general Gronwall's weighted inequalities; 
see Henry's famous book [18, Ch. 7]. 

Lemma 2.1. Let, for all t > 0, 

oo 

(2.15) JJ{u,r<Cse-K 

t n 

Then there exists a constant Cg > such that 

(2.16) \\u{t) - u{t)\\ < Cge"^ for all <t<T < oo. 
Proof. Firstly, if |t — r| < 1 we apply the Holder inequality to get that 

Mt) - u{t) r = j{]u,y<{r-t)j{ ]\u,f) < Cs e-*. 
n t Q t 



If |t — t| > 1, then we perform a uniform partition of the interval [t, r] into K parts of 
length 1 with a reminder, and apply the Holder inequality in each subinterval to get 

i=0 n t+i n t+i+K 

(2.17) 

□ 

Step 7: exponential estimate implies convergence to ip. We continue to describe 
the evolution of u{t) on a large finite interval + T] with A; ^ 1, on which, by 

assumption, the orbit enjoys all the estimates following from the optimal approximation 
on the stationary subset S^- So for any t E [t^, + T], there holds: 

(a) By (imi) and (I2A0|) . 

oo 

/ / {utY = B — V{u{t)) is exponentially small; 

t n 

(b) By Lemma [2. H — u{t)\\ is exponentially small; 

(c) Therefore, u{t) ~ u{tk) ^ for all t > tk, completing the proof in the case (i). 
The case (ii) is the same once we have fixed a uniformly bounded orbit. 

(iii) See (b) above. □ 

3. Discussion: resonance branching and exponential convergence 

3.1. Example: non-isolated equilibria and resonance branching. We take the 
equation with the cubic analytic nonlinearity in the unit ball in 

(3.1) AV^ + i/;3 = in Si = {|x| < 1} C M^, ijjldB, = 0. 

Then p = 3 is in the subcritical Sobolev range if 3 < f^, i.e., we need N < 4. Taking 
other p < ps provides us with similar examples for any N > 2. 

Firstly, we consider (13.11) in the half of the ball, i7+ = Bi (1 {xi > 0}. There exists a 
unique strictly positive classical solution -0+ G Hq{Q^) of (13.11) in constructed by the 
standard variational technique, [2| [29 j. 

Next, since equation (13.11) is invariant under the reflection xi ^ —xi, setting -0 = i/j^ 
in f2+ and ip = —"0+ in Q_ = Bi\ 0+ yields a non-radial i/g-solution of (13.11) in Bi. This 
extended function is a weak solution of (13. ip in the sense that it satisfies the corresponding 
integral identity and hence is a classical solutions by the theory of uniformly elliptic 
equations. Any smooth invariant orthogonal transformation in leaving the Laplacian 
and the Dirichlet boundary condition invariant produces continuous families of solutions 
consisting of non-isolated points. By the classical branching theory [20], in this case, 
belongs the spectrum of the corresponding linearized operator V"{i(j) = A + Si/j'^I. 
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In particular, we fix — 1 angels of rotations {6*1, ^jv„i} of the polar coordinate 
system in to get at least — 1 linearly independent elements of the kernel of A + 3-0^/ 
given by 

(/Pfc = ^^, A;= l,..,iV-l. 

Bearing in mind these rotations, the stationary subset 5'Ar_i is expected to be precisely 
(A^— l)-dimensional, and this k = N — 1 coincides with the kernel dimension. Unfortu- 
nately, we cannot prove the exact equality, and present later on another non-local model, 
for which Hypothesis (Hr) is guaranteed. For the problem on the plane (A^ = 2), the 
corresponding eigenfunction of V"{il)) with A = is given by rotation by the single angle 
e = 9i, so 

(3.2) ipi = ^^P, 

and the kernel seems to be one-dimensional (we do not prove this either). 

Using similar reflections, one can construct in the unit ball Bi C the nonlinear 
eigenf unctions ipilx) satisfying (13. ip with arbitrarily number / > 1 of connected positivity 
and negativity components. By orthogonal invariant transformation, each one generates 
continuous families of other stationary solutions. 

3.2. On exponential convergence. We first illustrate the reason for the exponential 
convergence by using the above example in M^, where the kernel is one-dimensional with 
the basis function ( 13. 2p . 

Orthogonality conditions for branching. The branching occurs in (13. ip from solution 
ifjQ = -0, SO by the classical branching theory [33], we look for a smooth curve of solutions 
in the form 

(3.3) i>^l = i>o + li-ipi + ^^^tp2 + + ^J'^i^4: + ••• , 

where, up to scaling, /i can be attributed to the angle of rotation. Substituting (13. 3p into 
(13. ip and equating similar terms yields the system for expansion coefficients 



(3.4) 



fAV'o + V'o = 0, 
A^i + 3V'o'^i = 0, 
AV^2 + ^Voi^2 = -3^0^?, 

A'04 + SipQip^ = -3V'o^2 - 3V'iV'2 - Gipoi'i'ips, 



The second equations says that ipi is from the non-empty kernel and hence is given by 
(13. 2p . The third equation then yields the first orthogonality condition 

(3.5) I^oi^f = 0, 

n 

which is necessary for branching to occur. Determining ?/'2 up to a constant (32, 



= ^2 + 
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where "02 is a solutions of the inhomogeneous equation, and substituting into the fourth 
one, for its solvabihty we have to have that 

Jijt + 6j Vo^?^2 + 6/32 / V-o^? = 0. 
n Q n 

Since the last term containing (32 vanishes by (13. 5p . we arrive at the second orthogonality 
condition 

(3.6) I^t + Qj ^o^!^2 = 0. 

n n 

Next choosing similarly 

V's = V's + Psipi, 
we obtain from the fifth equation that 

2 / i^0^f^3 + J V^O^l^i + / ^!^2 + I /?2 / = 0. 

n n n n 

This determines /?2 and so on. Thus, in general, this kind of branching demands more 
than one orthogonality condition (and, of course, existence of a non-trivial kernel). 

On a centre manifold link. Let us now discuss how these branching conditions affect 
the evolution properties of the corresponding parabolic equation 

(3.7) ut = Au + u^ in5ixM+. 

In this analytic case, there exists a local centre C°° manifold Wy^^{ipo)j see general theory 
in [26], [211 § 9]) [32] for applications to reaction-diffusion equations and [15], where the 
convergence problem is studied. 

Since Wi^^^ipo) is tangent to the corresponding centre subspace i?^(0) = Spanl-^i}, 
looking for the centre manifold behaviour with solutions of the form (see references and 
accurate estimates for the non-analytic flow in the next section) 

u{x,t) = iJjq{x) + ai{t)(f)i{x) + ... for t > 1, 

and projecting the PDE (13.71) on E'^{0), we obtain a perturbed ODE 

(3.8) a[ = -fal + ..., 
where, as one can see, the constant 7 vanishes, 

(3.9) 7 = /^oV'? = 0, 

by the first orthogonality branching condition (13.51) . Therefore the higher-order terms 
should be taken into account, but these seem also do not help to detect a suitable non- 
trivial evolution on the centre manifold. 

Indeed one can see that this one-dimensional branching manifold 5*1 is a centre manifold 
for the parabolic problem. Under the above assumptions, denoting by {(pk} the complete 
subset of eigenf unctions of V"{iPq) = A + Sip^I, we may assume the following expansion 
(cf. (O) 

Si = {i>^, = ipo + ai(/i)0i + a2(/i)02 + •••}, 
11 



where 0^(0) = for any k. Therefore, the evolution on this centre manifold with 



u{x, t) = i/jQ^x) + ak{t)(f)k{x) 
(k) 

is governed by the trivial dynamical system 

4 = 0, fc = l,2,..., 

which we have observed in (I3.8p . (I3.9p . Notice that a centre manifold consisting of pure 
equilibria is not an exceptional situation in parabolic problems; cf. [8], where such an 
invariant exponentially stable centre manifold has been detected in the free-boundary 
MuUins-Sekerka model (a kind of Hele-Shaw flow). 

Thus in this well as and in other resonance cases, there exists a local centre 

manifold for the parabolic problem precisely coinciding with the stationary subset S^, so 
no evolution of orbits on W{^^{iIjo) can be observed. Then the exponential convergence in 
the Theorem ll.il (iii) can be associated with the fact that, in the orthogonal complement 
of the kernel, the behaviour is purely exponential and corresponds to the evolution on the 
orthogonal stable manifold constructed at a different equilibrium = ipf_i(t) ~ ^o! cf. Step 
1 of the proof in Section [21 In the invariant manifold theory, this is usually expressed by 
the fact that, for nonlinearities and good linearized sectorial operators with discrete 
spectrum and finite-dimensional unstable and centre subspaces, etc., the centre manifold 
is exponentially stable; see e.g. [HI Prop. 9.2.3] and [32] . 

4. An explicit example with non-local nonlinearity 

Let Q he a bounded smooth domain in and let {A^} and {ipk} be the eigenvalues 
of A, where each one repeated as many times as its multiplicity > 1? and the cor- 
responding complete, closed subset of eigenfunctions that are orthonormal in L^(f2). The 
case of the unit ball, Q = Bi, is classical. Here, for A, all the eigenvalues, eigenfunctions, 
multiplicities, etc., are described by the Laplace-Beltrami operator A^. on the unit sphere 
S*^"^ = dBi in M^, which in the polar coordinates (r, a) is given by 

(4.1) A = A, + ^A<„ 

Act is a regular operator with the discrete spectrum in L^(S'^~^) (again each one repeated 
as many times as its multiplicity) 

(4.2) (^i-K) = {i^k = k{k + N ~2), k> 0} 

and an orthonormal, complete, closed subset {Vfc(cr)} of eigenfunctions, which are homo- 
geneous harmonic k-th order polynomials restricted to S^~^. 

Fix an / > 2 and consider the PDE with a non-local cubic nonlinear term 

(4.3) ut = Au + Xiu- {Ju^)u 

with the same Dirichlet boundary conditions and bounded initial data Uq. Such non-local 
parabolic models can provide us with further examples of delicate asymptotics obtained 
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via explicit computations; see [29]. Firstly, it is easy to solve the stationary equation 
(4.4) A^ + Az^- (/V'')^ = 0. 

Namely, there exist non-trivial equilibria of the form 



(4.5) ■0(x) = ± a/Az — Xj ifjj (x) provided that A/ > Xj, 

and other obviously constructed linear combinations of such functions corresponding to 
the same eigenvalues. 

4.1. Result on convergence. Secondly, studying the parabolic equation (14.31) and look- 
ing for solutions in the form of the eigenfunction expansion 

U{x,t) = Y.O'k{t)'lpk{x), 

(k) 

we obtain the following dynamical system for the expansion coefficients: 
a'), = {Xi — Xk — \a\'^)ak, fc = l,2,..., where \a\'^ = J2'^k- 

Integrating yields ak{t) = ak{0)e^ '''^e o . Calculating the sum |ap, we derive 

-}\a\^(s)ds 

|a|2 = ^(t)e , with ^(t) = ^e2(^'-^^)*ai(0). 

Setting Z = e-Zo yields a simple ODE, Z' = -2^(t)Z2, so, on integration, 

af(0)+ E c(^i-^'")*a2„(0) 

where for simplicity af{0) denotes the sum ^ ctfc(O). Finally we arrive at the following 

Afc=A; 

expressions for the expansion coefficients: 

(4.6) ak{t) = ak{0) 



This reveals three different cases of solutions with exponential and algebraic rate of con- 
vergence to equilibria, depending on the dominant term in the long square root in the 
denominator. Namely, we take initial function 

^o(a;) = J2 afc(O)'0fc(x) 
such that there exists j > 1, for which 

(4.7) ai(0) = ... = aj_i(0) = and aj{0) ^ 0. 

It is easy to derive from (14.61) the following result. 
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Proposition 4.1. Let ( |^. 7p hold. Then, as t —>■ oo, 

(4.8) (i) M(x,t) = 0(e-(^^-^')*) ^ 0, ifXj>Xi; 

(4.9) (ii) M(2;,t) = O(^) ^ 0, if Xj = Xi; and 

(4.10) (iii) u{x,t) = tlj{x) + 0{te-^^^'-^^^'), if Xj < Xi, 
where ip ^ is an equilibrium. 

It is important that the case (iii), where branching of equihbria is available by (14. 5p . 
precisely demonstrates that the convergence to non-trivial stationary solutions is expo- 
nential; cf. Theorem II. l( iii). It is worth mentioning that the rate of convergence in (14.101) 
is not purely exponential and, in general, contains a lower-order algebraic factor t that 
is induced by the C'(t)-term in the square root in (14. 6p . Different In t-perturbations can 
occur in problems with usual (local) nonlinearities; see the Remark in Section O 

The only case in Proposition 14. H where the convergence is not exponentially fast, is 
(ii), in which, by (14.51) . the only equilibrium is trivial, ■0 = 0, and is isolated. 

4.2. Hypothesis (Hr) is valid. We now present a rigorous evidence that Hypothesis 
(Hr) is valid for such operators. Consider the elliptic equation (14. 4p . Fix a j > I, let 
Kj = 1 + k > 2he the multiplicity of Xj and let ...,ipj,k+i} be the orthonormal subset 
of eigenfunctions of A corresponding to Xj. We fix an equilibrium 

fc+i 

(4.11) ^{X)= Ec^iJJA^), 

i=l 

where, on substitution into (14.40 . the coefficients satisfy (cf. (14.50 ) 

(4.12) Ecl = Xi-X,. 

Stationary subset 5*^. Obviously the subset containing given equilibrium (14. lip is 
described by 

fe+i 

(4.13) i; = 'ij){x) = Ciipj^i, where ^c,^ = A; - Xj. 

i=l 

Therefore 

(4.14) dim Sk = k = Kj -1>1. 
Kernel of V"{i>). It follows that 

(4.15) V"{ip)v = Av + Xiv-2i){J'4jv) -v{J . 
Substituting (14. lip yields the following equation for the kernel: 

(4.16) V"{'4j)v = Av + XjV - 2tp{J ipv) = 0. 
Finally taking v in the form 

fc+i 

i=l 
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and substituting into (I4.16p gives a single condition on k+1 expansion coefficients 

fc+i 

E = 0, 

i=l 

SO that the kernel is precisely fc-dimensional. 

By (14.141) this completes the analysis and proves that (Hr) is always valid for such 
non-local elliptic operators. Without any changes this example extends to such operators 
of arbitrary order 

V'{t/j) = -(-A)™?/) + XiiIj- {J with any m > 1. 

5. Rate of convergence can be arbitrarily slow: ^ = 

Here we present some estimates showing that the rate of convergence in parabolic 
problems with specially designed nonlinearities can be arbitrarily slow. 

5.1. The original semilinear parabolic equation. Let f2 be a bounded smooth do- 
main in such that Ai = —1 is the ffist simple eigenvalue of A in L^(f2) with domain 
Hq{Q) and the normalized eigenfunction 0i > in Q. Consider the following parabolic 
equation: 

(5.1) Ut = V'{u) = Au + u-e-^/''^ in x R+, u = ondQxR+, 

where we set /(O) = by continuity. Here e~^/"^ is the standard function, which 
is not analytic at m = 0. In view of the perfect spectral properties of A, we have that 
there exists a local invariant, C°°, one-dimensional centre manifold ^2^^,(0) of operator 
V'{u), which is tangent to the eigenspace E^{0) = Span{0i} of A-|-/; see [211 Thm. 9.2.2] 
and [32]. Looking for the corresponding center manifold behaviour, we decompose the 
solution in the form 

(5.2) u{x,t) = ai{t)(f)i{x) + U2{x,t) for t ^ 1, 

where U2{x,t) = ^ ak{t)(j)k{x) ± E'^^O), U2{-,t) = o{ai(t)) as t ^ oo and we assume that 

k>2 

ai{t) > for t ^ 1, e.g., we take positive solutions. Projecting the PDE onto £"^(0) yields 

(5.3) a[ = -(e-i/"',0i) = -(e-VK<^?+-), 

for t ^ 1, where (■, ■) denoted the inner product in L^(f2). Denoting pi = max(/)i(x) > 
and P2 = / 01 > yields the ordinary differential inequality 

(5.4) a[ > -p2e-i/^^i"?, 

so that ai{t) > ai(t), where ai solves the corresponding ODE with the equality sign in 
(15. 4p . Finally, we obtain the following estimate: 

(5.5) «iW>5^ ^ort»l 

for the stabilization on the centre manifold. Obviously, this corresponds to ^ = in 
(11.141) . so that the rational algebraic modulus | ■ |^~^ in (I1.12p cannot be applied to such 
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C°° nonlinearities. One can introduce an appropriate modulus, for which the rate (15 .51) 
is acceptable. For instance, for a slightly modified function 

we have the generalized gradient inequality inequality at -0 = for u ^ (actually, it is 
equality) 

u{V{u)) < \V'{u)\ with modulus u{s) = 2\s\ \ ln(2s)|^/^ 

which is "almost" linear, strictly concave function as s ^ 0. 

It is easy to present other C°°, non-analytic nonlinearities generating arbitrarily slow 
rate of convergence. Recall that, according to the results in Section [31 any such slow 
convergence corresponds to isolated equilibria. 

Remark: a cubic nonlinearity. In the presence of non-trivial kernels, a full asymptotic 
expansion of solutions can be a difficult problem even in the analytic cubic case 

(5.6) Ut = Au + u-u^. 

It was shown in [1] that, under the same kernel assumption 

(5.7) dimker (A + /) = 1, 

there exist solutions with the following logarithmically perturbed decay as t — > cx): 

n k 

uix,t) = E t-^-' E ^kj{x){\nty + o{t-2-^), 

k=0 j=0 

for some integer n > 0, where ipkj are solutions of well-posed elliptic problems. This 
corresponds to a special case of centre manifold behaviour, where the leading term of 
convergence to = is of order 0(t^2) that does not contain the Int-factor. It is 
easy to see why in this case branching of equilibria from zero is not possible. Indeed, 
if this occurs, then in view of (15. 7p . by branching theory [33], the equilibrium branch 
should have the representation ip = ^ipQ + where ^ is the branching parameter and 
(fo is the eigenfunction of A with A = — 1. Substituting this into the stationary equation 
Alp + -0 — -0^ = and multiplying by Lpo yields J^ip^ = 0, a. contradiction, so -0 = is the 
isolated equilibrium. 

5.2. A non-local semilinear parabolic equation. As usual the computations are sim- 
plified for semilinear non-local parabolic fiows like 

(5.8) Ut = Au + fiu — g'(^J u'^)u, 

n 

where g'{s) is a given C°°, non-analytic function. The potential here is 

(5.9) V{u{t)) = -I / \Du{t)\' + llJu'-g{Ju'). 

Assuming again that /i = Ai is the first eigenvalue of A in L'^{Q), we obtain that the 
centre manifold behaviour (15. 2p is described by the ODE 

(5.10) a[ = -g{a\ + ...) for t > 1. 
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Choosing non-analytic C°° functions g, e.g., 

g{s) = s^/^e-i/^ for s > 0, 

integrating (15.101) asymptotically yields a non-algebraic decay 

diit) ~ „ as t ^ oo. 

Notice that, besides slow decay behaviour, the present "less nonlinear" model (15. 8p is 
suitable for revealing refined evolution properties of orbits describing stabilization phe- 
nomena. 

6. On applications to smooth gradient systems in Hilbert spaces 

Such an extension is straightforward, since, in the proof of Theorem 11.11 in Section [21 
we have minimally used the specific properties of the second-order elliptic and parabolic 
equations under consideration. 

Denoting by (•, ■) and || ■ || the inner product and the induced norm in a separable 
Hilbert space if, we consider a smooth gradient flow in if, 

(6.1) ut = V'{u) for t > 0, 

where the operator V with a dense, compactly embedded domain = D{V') G H is the 
Frechet derivative of a C^-functional V : M. We assume that V"{u) : — >■ if is a 

Fredholm operator of index zero admitting a suitable self-adjoint extension with discrete 
spectrum, compact resolvent and a subset of eigenfunctions, which is complete and closed 
in H. We impose the necessary condition of coercivity of V to guarantee existence of 
global orbits {u{t)} for arbitrary initial data uq G H^, which are sufficiently regular to 
satisfy the gradient identity (cf. (I2.10p ) 

(6.2) ^^V{u{t)) = \\utr>0. 

We impose the same main Hypothesis (Hr) and will next check, using the scheme of 
the proof from Section [21 which conditions we need to guarantee the result (i) and, hence, 
(iii) in Theorem 11.11 

Steps 1—3. The construction remains the same and (12.21) . (12. 5p and (12.71) are guaranteed 
by assumed good spectral properties of the linearized self-adjoint operator V'^ipiJ.) ^oi any 
tfj^ E Sk uniformly in /i ^ 0. 

Step 4. Estimate (12.91) is valid for sufficiently smooth operator V in if^. 

Step 5. Here we need estimate (]2.12p . which is indeed Lagrange's formula for the smooth 
functional V in H^. 

Step 6. Lemma [2.11 remains valid in the topology of H. 
Step 7 remains unchanged. 

Hence Theorem ll.il (i). (iii) is true for arbitrary smooth gradient flows under the presence 
of resonance branching. It follows that the convergence (evolution completeness) result 
holds for the 2mth-order PDEs (I1.16P and for other smooth gradient parabolic flows 
governed by higher-order operators. 
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7. Non-autonomous perturbations of gradient systems 

More modifications of the approacli are necessary to prove tlie convergence result in 
Tlieorem 11.11 (i) for non-autonomous perturbations of (16.11) . The main features of our 
analysis are illustrated by the following example: 

(7.1) ut = V'{u) + h{t)W'{u) for t > 0, 

where W{u) : R is a C^-functional. Without loss of generality and for simplification 

of future manipulations, we assume that the decay rate of perturbation h E satisfies 

(7.2) h{t) as t ^ oo, and h{t) > 0, h'{t) < for t > 0, 

so that, in particular, h' G L^(]R+). Otherwise, for non-monotone and changing sign 
rates of perturbations, we can use estimates from above and below with functions h± (t) 
satisfying necessary assumptions. 

Actually, we will need a more restrictive condition on the decay rate, 

(7.3) VheL\R+). 

It is important to deal with non-exponentially small perturbations such as 

(7.4) h{t) = jY+t)^ '^i^h a > 2, 

for which (17. 3p holds. For the future purpose of integration of ordinary differential in- 
equalities, we characterize this class of such slow decaying functions as follows: for any 
constant (3 > 0, 

t 

(7.5) //i(s)e^Ms = i/i(t)e^*(l + o(l)) as t ^ cx), 



or, equivalently, by L'Hospital's rule, 

(7.6) ^ °~ as t ^ oo. 

Since (17.11) is not a gradient system in the sense that a monotone Lyapunov function, 
in general, does not exist, we need first to prove that the crucial characterization (II. 6p of 
the w-limit set remains valid, so we should begin with 

Step 0: u{uo) C $. Multiplying (17.11) by Ut in H, instead of (16.21) . we obtain the identity 

(7.7) A v{u{t)) + ^ mW{u)) - h'{t)W{u) = \\u,r> 0. 
Integrating over [t, oo) and using that h' G L^(M+) yields the convergence 

oo 

(7.8) J\\utis)fds<oo. 

t 

Hence, given a sequence {tk} oo such that u{tk) — > "0 G uj{uo), we obtain by the Holder 
inequality that, for arbitrarily large fixed t > 0, 

CX) 

(7.9) \\u{tk + t) - u{tk)f < t J \\ut\\^ ds ^ as/c^cx). 

18 



Passing to the hmit as t = tk + 1 oo in (17.11) . we then conclude that u{tk + 1) converges 
in i^^ to a solution u{t) of the limit autonomous equation (16.11) . which is independent of 
t, so it is an equilibrium ip E ^. 

We now return to seven steps of the proof in Section [21 

Steps 1-3 are unchanged. 

Step 4. In view of the extra term in (17. ip . instead of (12. 8p . we will have 

(7.10) g = ut + giiu) + h{t)W'iu). 
Therefore, by (17.20 . instead of (12. 9p . we obtain holds 

(7.11) Mt)-^^^t)\\H2<C,i\\ut\\+h{t)). 

Step 5. It follows from (17.70 that the convergence (12.110 takes place but not necessarily 
from below, which is not important. Then, instead of (I2.13p . we have 

(7.12) f^[B- V{um = -Wutf + {h{t)W'{u),Ut) < -Cr[B - V{u{t))] + Cjh{t). 

Integrating this inequality by using the slow decay hypotheses (17.50 . (17.60 . we conclude 
that, instead of the exponential estimate (I2.14p . on large bounded intervals t E [tk,tk + T] 
a.s k —>■ oo, with T ^ 1 and t 3> 1, 

(7.13) B -V{u{t)) <2h{tk + t). 

Step 6 needs a major revision, since Zelenyak's Lemma 12.11 applies only to exponential 
decay estimates. 

Lemma 7.1. Assume ( |7.^P and l \7.3^ hold and, for all t > 0, 

oo 

(7.14) J\\ut{s)rds<h{t). 

t 

Then 

oo 

(7.15) \\u{t) ~ u{t)\\ < J y^h(^ds for all <t <r <oo. 
Proof. The result is obvious if |t — < 1, where, by the Holder inequality. 



\u 



(t) - u{r)\\ < J J \\ut{sWds <VT^t^h{t). 



For large intervals t — t > 1, we perform a partition with a sequence of time-steps 
{Ai,i = 0, 1, K}, so that, instead of (I2.17p . we obtain 



K oo K 

(7.16) ||n(t)-n(r)|| <E v^W/lkt(5)Pds = E 

i=0 V ti i=0 



where tj+i = + Aj. Setting Aj = 1 for i = 0,1, K — 1 and A/< < 1, by the integral 
test of convergence of series, we obtain (I7.15p . □ 
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Step 7 remains the same, where we replace "exponentially small" in (a) and (b) by 
"C»(v/M^)-smair. 
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Comment . Almost two years after publishing the present paper irl^ 

Sbornik: Math., 198:6 (2007), 817-838; see [H], 

in January 2009, the first author found that similar Zelenyak's ideas and techniques were applied 
in \17\ § 5] to convergence in ID wave equations. Actually, Zelenyak's Lemma [2. II were introduced 
in [T7] earlier, in Section 2, and was a key ingredient of the proof of Theorem 1.1 on exponential 
convergence in a Hilbert space setting (similar to our Theorem 1.1 in an analogous Hilbert space 
framework explained in Section [6|) . It is interesting and truly remarkable that such a growing 
interest to Zelenyak's fundamental idea on proving exponential convergence from the 1960s |34j . 
in the 21st century, occurred approximately in the same year, about 2007, according to the 
publication dates. 
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